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Abstract. We show that the formalism commonly used to implement Bose- 
Einstein correlations in Monte-Carlo simulations can lead to values of the two- 
particle correlator significantly smaller than unity, in the case of sources with 
, strong position-momentum correlations. This is more pronounced when the 

' phase space of the emitted particles is strongly reduced by experimental accep- 

ts) , tance or kinematic analysis selections. It is inconsistent with general principles 

' according to which the Bose-Einstein correlator is larger than unity. This in- 

consistency seems to be rooted in the fact that quantum mechanical localization 
' properties are not taken into account properly. 



1 Introduction 

Hanbury-Brown Twiss interferometry has been used |2| in both high energy and nuclear 
physics to determine the space-time dimensions of the emitting source created during nuclear 
collisions by using the effect of the interference pattern between two identical produced bosons 

iSI- 

The source parameters derived from fits to the correlation function are difficult to interpret 
, directly as real geometric quantities, being sensitive to the transverse and longitudinal dynamical 

' expansion of the system B f?!, |^ , ||, |o[ 0, which result in a momentum dependence of the 

extracted source radii |]l2| , [l^, |14[ to long lived resonance decays @ |7[ IH, @ , the Coulomb 
interaction and final state rescattering (22j. There has been recent experimental evidence 

for flow effects and one possible implication of flow is that distant points of emission in the 

source volume cannot emit particles with closely differing momenta, and thus do not contribute 
to the small relative pair momentum region |Q. It is also anticipated that strong absorption 
must exist in the case of large stopping power; a particle originating at the side of the source 
opposite to the direction of its momentum cannot easily propagate through the source to be seen 
by the detector, and therefore only a limited region of the source will be seen, noted already in 
AGS studies |25j. It is thus interesting to try to further probe the relationship between source 
geometry, dynamical expansion, and kinematical regimes viewed by the measuring apparatus. 

To study the effect of position-momentum correlations on the shape of the correlation func- 
tions, a simple Monte Carlo phase-space model controlled by a few macroscopic parameters was 
developed. As is the case with more detailed and sophisticated microscopic event generators, 
there is no Bose-Einstein symmetrization effect included from first principles j2^, . The 

Bose-Einstein correlations were then added to the initial distributions by including the sym- 
metrization in the form of a weight calculated for each pair of identical particles, a procedure 
found extensively in the literature j2^, |2^, 29 



In the present work we show that there are limitations to this formalism, as it is an ap- 
proximation and several assumptions are implicit in its derivation. This will be demonstrated 
on the basis of the Monte-Carlo calculations. We briefly check that our model produces results 
consistent with theoretical predictions as well as with other simulation studies in case of similar 
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phase space distributions. In our further studies of source models with position momentum 
dependence, we then find that the current praxis of including Bose-Einstein correlations is in- 
consistent with general principles. Hence we restrict our discussion to pointing to a possible 
origin of this problem. 

This paper proceeds as follows: Section 2 describes the Monte-Carlo model and the inclusion 
of the Bose-Einstein effect, Section 3 presents the model results and in Section 4, we turn to the 
discussion of the observed problems. 

2 The Model 

Our model produces events from an analytically given phase space distribution by a Monte Carlo 
technique. The model provides the phase space distribution of particles at the points [x^ ,p^) of 
their last interaction, with no assumptions about the dynamical evolution of the collision. The 
collision region is described in terms of a few macroscopic parameters defining the spatiotemporal 
extension of the source, such as the source shape and size, and the dynamical features of the 
system, such as temperature and collective flow. All the particles are assumed to be pions, and 
resonances are not included. 

As for all existing event generators, the obtained phase space distribution does not contain 
Bose-Einstein correlations. To include the latter, we applied the following prescription, widely 
found and used in the literature |2^, |3^. Each identical pion pair emitted from the points 
(ri, ti) and (r^, is weighted with the Born probability of a symmetrized two-pion plane wave, 

I * 1 + cos\{r-{ - r2) ■ (pi - V2) - {h - t2){Ei - E2)] (1) 

which reduces to 

I * |2= 1 + cos[(7-l - r1) • (pi - p'2)] (2) 

since the pions at freezeout are on-shell by definition and the source emission is assumed to be 
instantaneous (there is no time evolution) . The correlator defined via Eq. (||) is a function of the 
three-dimensional relative momentum component q. The correlator C(q) is then obtained for 
each bin as the sum over pion pairs weighted by | 4* ^ and normalized to the sum of unweighted 
pion pairs, cf. The Coulomb interaction between the pions is not simulated. 

Here, we streamline our presentation by restricting it to the results of the one-dimensional 
fits of C((?3) in terms of the 3-momentum difference 

Ciqs) - 1 + Xexp[-qlRl] qs =| p{ - p'2 \ ■ (3) 

All calculations are done in the longitudinal comoving system LCMS frame, which is defined 
such that the longitudinal component of the average pair momentum vanishes. 

3 Model Results 

In contrast to full event generators like Venus, RQMD, ARC, etc., which try to incorporate all the 
physics expected to be present in a heavy ion collision, the purpose of our model was to isolate and 
study one important effect: the geometrical and dynamical interpretation of HBT parameters 
in the presence of radial flow and realistic experimental acceptance. The model's simplicity 
allowed the well controlled study of a wide set of different flow and acceptance conditions. Thus 
in the course of this study, we have found, inadequacies in the common practice (Eq. (^) 
of including Bose-Einstein correlations which become particularly apparent for sources with 
strong position-momentum dependence where certain kinematical selections are imposed. Our 
presentation will first illustrate this effect in an instructive way for the extreme case of complete 
position-momentum correlations in the source. Next, we will summarize the results of flow and 
acceptance effects on the HBT parameters for realistic phase space distributions and acceptance 
criteria. 
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3.1 Source with Complete Position-Momentum Correlations 

Here, we calculate the correlation function for a linear source in the beam (z) direction with 
neither transverse spatial extension nor transverse momentum dependence. The longitudinal 
momentum is chosen to be completely due to flow, 

piz)^Dz, (4) 

where p has units of GcV/c, D being a constant. This distribution represents a source expanding 
in the z-direction for which the argument in the cosine of Eq. (^) reduces to (rl — t^) ■ (pi —pi) = 
Q\/D^. For D — 0.02, this leads to the correlation function shown in Figure |l|. The source 
in Eq. (^) has a total position-momentum correlation and in this case, as can be checked 
analytically |^ , the correlator obtained for Eq. (|^) with the cos-prescription oscillates between 
2 and 0. Introducing a Gaussian spatial smearing of the emission points in Eq. as might 
be motivated by the picture of a limited quantum mechanical localization of particles, one sees 
that the oscillations of the correlator decrease. Still, the correlator can drop below unity, and 
this is in strong contradiction to calculations from first principles p6| which ascertain that for 
arbitrary sources the correlator is always larger than unity. We next investigate in how far this 
behaviour persists for more realistic phase space distributions. 



3.2 Position-Momentum Dependence in a Realistic Model 

To study more realistic scenarios, we generate pions according to macroscopic model parameters 
that correspond to an instantaneous Gaussian source with realistic momentum distributions. No 
attempt is made either to reproduce any data or to create a fully realistic Monte-Carlo event 
generator. Accordingly the spatiotemporal part of the source is modelled by 

G = const. X exp [-{x^ + + z^)/R^] ■ (5) 

The transverse momentum dependence and the rapidity dependence chosen for the case of no 
flow (i.e., no position momentum dependence in the source) are 

dN/dpT — ApT ey.^ {~pt / B) , 

y = c%/-21na(cos(27r6)) . (6) 

As an input parameter, we use a Gaussian radius R — 6 fm, motivated by the hard sphere radius 
of the ^^^Ph incoming projectile. The input for the momentum distributions is chosen according 
to the measured transverse momentum (j)t) and rapidity {y) distributions in the CERN 158 
GeV/n Pb-|-Pb data The inverse slope parameter of the pT-distribution is taken to be 

B = 200 MeV and A is an arbitrary normalization constant. The rapidity {y = \ log ^'^^ ) 
dependence is specified by random numbers a and b which arc uniformly distributed in the 
interval (0,1), c being a constant. 

To incorporate radial flow in the model, we modify the phase space distribution (Eqs. 
by introducing a radial flow (3{r) of the emission points, 

/3(r) = 1 - e-^'f , (7) 

where / is an adjustable parameter. For different flow strengths /, the radial dependence is 
shown in Fig. ^. Superimposed is the mean value of the radial velocity 

A = #|^- (8) 

extracted from simulated Pb+Pb Venus events, version 4.12 |2^. Here E denotes the total 
energy of the pion and r its radial distance from the source center. One sees that a choice of 
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/ = 9 fm fits the Venus data very well. This value of the flow strength, / = 9, was used for 
simulating a "realistic" flow. No attempt was made to reproduce any other Venus distribution. 
As demonstrated in , the flow dependence of the observables is mainly due to the size of the 
flow, while its functional shape plays a somewhat secondary role. The order of magnitude of the 
flow velocity extracted here 0.35c) can be compared to the flow parameter rjf in [p^ . 

We have generated events for different values of the flow strength / and the result was 
verified to be independent of statistics. Each Monte-Carlo event contained typically 100 identical 
pions. For the relative pair momentum differences, 5 MeV bins were used and the results were 
checked to be insensitive to bin sizes in the range of 5-20 MeV. The original aim was the 
study of the flow and acceptance dependence of HBT radius parameters. Indeed, our simple 
model shows reasonable physical properties. Especially the flow dependence of the 1- and 3- 
dimensional kr integrated HBT radius parameters is in qualitative agreement with that obtained 
in other model studies, |^5[ |2^: The HBT radii decrease with increasing ffow strength /, since 
the effective emission region ("region of homogeneity" [p^) for pion pairs with small relative 
momenta decreases for increasing /. Also, we have considered the so far little studied effect of 
the detector acceptance on the HBT radius parameters. In Figure || (left), we show for the case 
of a realistic flow strength / = 9 fm the freezeout positions of all generated pions and of those 
satisfying the {px , y) acceptance criteria of a typical magnetic spectrometer [|^ . The emission 
region of the detected pions is clearly smaller than the total emission region. The corresponding 
HBT radius parameter R3 is presented in Figure ^ (right) as a function of the flow strength / 
for the cases with a realistic magnetic spectrometer acceptance and without. In both cases, the 
HBT radii decrease with increasing flow and the acceptance dependence is very small. ^From 
all this we conclude here only that our model is not too oversimplified and reproduces essential 
features obtained in more complete simulations. 

We now consider a simple modification of our model, the introduction of an absorption cut. 
To this aim, we impose on the Monte Carlo output, defined by the Eqs. (||-|1), a kinematical cut 
which effectively strengthens the position-momentum correlation in the source: pions are only 
emitted, if their momentum vector is in an angle less than 45 degrees around the direction of 
their position vector, i.e., each pion has to move out of the source in a 45 degree cone. Under 
these conditions, the correlation functions in 93 displayed in Figure |4|have been obtained, where, 
additionally, cuts in single particle transverse momentum of 0.1 GeV/c - 0.3 GeV/c have been 
applied. 



We emphasize that in contrast to the pathological source discussed in subsection 3.1, this 



modified source shows a rather reasonable phase space distribution which is difficult to reject a 
priori. Still, the correlator obtained from the cos-prescription (Eq. |l|) again drops below unity. 
Especially, when px is restricted to small values, the dip in the correlation function around 0.08 
GeV/c figures more prominently. Such a characteristic dip can be produced in C^q^) given at 
least a px cut, and either radial flow or the angle cut. We thus conclude that it arises when 
imposing an acceptance cut on the pions emitted from the Gaussian model with strong position- 
momentum correlation. It is worth noting that the half- width of the correlation function remains 
nearly identical in all four cases in Figure 0. 



4 Discussion 

Existing event generators do not propagate (anti)-symmetrized wave functions and hence face a 
conceptual difficulty in incorporating the effect of Bose-Einstein correlations |Q . The current 
practice of modifying the weight of pion pairs by the Born probability \ ip p of symmetrized plane 
waves does not address this problem properly. Here, we have shown for the first time that this 
conceptual difficulty can have significant quantitative consequences. In fact, the dip observed 
in the correlator Fig. ^ is clearly unphysical and shows that the formalism leading to Eq. (||) 
breaks down. This becomes clear if one looks at the formalism used for analytical models of 
the emission function, where the ad hoc prescription (111) is not needed. In these, the evaluation 
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of the Bose-Einstein correlation function is based on the coherent state formahsm arising from 
quantum field theory, using p6[,p^ 

' jSxS{x,K+\q)^d^xS{x,K~\q) ^' 

with q = pi ~ P2 and K = (pi +p2)/2. Clearly, this correlation function cannot become smaller 
than unity. The approximate method using Eq. (^) is based on a semiclassical picture for a set 
of discrete space-time points. It produces effects that are not always consistent with Eq. 
and which can become non-negligible as seen in Figure ^. They are more pronounced in the 
presence of strong position-momentum dependence when the long-range characteristics of the 
argument {fl — r^) • {pi — P2) in @j play an important role. 

The discrepancy in methodology can be even more drastic; an analytical calculation based 
on Eq. (||) using the source given in Eq. (|^) does not show these oscillations while using 
Eq. with the same source yields Figure |l]. However, Eq. (|^) is at this time the only 
method available to build correlation functions from the space-time output of microscopic event 



generators |28 ,|]16[, , ||3l| , . This formalism has yielded reasonable results under the less 
severe conditions studied until now. The exact expression given in Eq. (^) however cannot be 
used in a Monte-Carlo simulation when the Monte-Carlo event generator does not provide the 
source density function of the mean momentum K. The probabilistic Monte-Carlo approach 
does not allow to deal with the quantum mechanics effects involved here. Without a solution 
to this problem, since the procedure based on Eq. (0) has been used by many groups in the 



recent literature |g8[ , Ijl^l , 1 22 , [g9| , it is now henceforth clearly important to improve the current 
Monte-Carlo Bose-Einstein formalism and to estimate the errors involved in the procedure using 
Eq. (I). 



5 Conclusions 

We have shown that the cos-prescription commonly used to include Bose-Einstein effects in the 
Monte-Carlo simulation can lead to results not consistent with first principles. This calls into 
question its quantitative and qualitative reliability, especially for the case that certain kinematical 
selections are applied when strong position-momentum correlations are inherent in the pion 
source. A deeper understanding of this simulation formalism is necessary in order to make 
more detailed analyses of dynamical issues and acceptance effects. Most remarkably, the cos- 
prescription (|l|) interprets both position r^, ti, and momenta pi returned by an event generator 
as sharp (classical) phase space coordinates. This clearly violates the Heisenberg uncertainty 
principle. On the other hand, a smearing of the emission points, motivated e.g. by the picture of 
a limited quantum mechanical particle localization, allows to remedy the unphysical dip in the 
correlator at least partly, see Section |3.1| . This observation may indicate in our opinion that the 
inconsistencies of the prescription (|l|) presented here are rooted in an incorrect treatment of the 
quantum mechanical particle localization. This points out the need for an advanced quantum 
mechanical Monte-Carlo event generator that can properly describe Bose-Einstein correlation 
functions. 
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Figure 1: For a linear source expanding in the z direction according to Eq. (j^), D = 0.02. Left: 
Momentum pz in the beam direction as a function of the z position of emitted pions. Right: 
The resulting two-pion correlation in using the formula in Eq. (0) . 




Figure 2: The radial flow velocity as a function of the radial distance from the center of a pure 
Gaussian source parametrized according to Eq. (||). The flow extracted from Venus (v. 4.12) is 
well represented by this parametrization with /=9 fm (filled triangles). The vertical line shows 
the true 6 fm Gaussian source size. 
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Figure 3: Left: Freezeout positions in the x direction for all pions and those satisfying a 
kinematical acceptance cut, where the initial pion momenta have been changed by adding radial 
flow extracted from Venus (eg. Eq. (||) with /=9 fm)). Right: The HBT radius R3 as a function 
of the flow velocity at a distance of r=6 fm from the source center for the radial flow profiles in 
Figure ||. 
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Figure 4: Pion correlation in (J3, with a horizon cut imposed such that the pion momentum 
vector must be within a 45 degree angle of the radial vector, for four different selections in single 
pion transverse momentum. 



